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Abstract 

This  paper  investigates  the  length  scale  of  ocean  surface  breaking  waves  m  the  spectral  range  of  intermediate  wavelength 
components  a  few  centimeters  to  a  few  meters  long*  The  spectral  properties  of  wave  breaking  are  examined  first  with  the 
dissipation  function  of  the  wave  action  density  conservation  equation.  The  analysis  reveals  a  strong  breaking  signature  in  wave 
components  between  0. 1 5  and  1 .5  rn  long  in  the  form  of  a  quasi-singular  behav  tor  of  the  dissipation  function  using  the  present 
formulation  of  the  wind-generation  and  breaking  dissipation  functions.  Independent  studies  of  more-direct  breaking  observations  of 
radar  tracking  of  sea  spikes  in  the  past  have  shown  close  correlation  between  sea  spikes  and  scattered  traveling  at  the  speed  of 
surface  waves  a  few  meters  long  and  much  shorter  than  the  dominant  w  avelength.  This  feature  of  sea-spike  properties  is  consistent 
with  the  breaking  signature  of  the  dissipation  function  in  similar  wavelengths.  The  intermediate-scale  waves  are  the  primary 
contributor  of  the  ocean  surface  mean-square  slope.  The  close  correlation  between  the  gas  transfer  rate  and  the  mean-square  slope 
has  been  demonstrated  repeatedly  A  better  understanding  of  the  wave  dynamics  of  intermediate-scale  waves  is  important  lor 
clarification  of  various  gas  transfer  mechanisms, 
v  2006  Elsevier  B.V.  All  rights  reserved. 

Ki'vwonts  Wave  breaking;  Surface  roughness,  I nlermedi ale-length  scale;  Dissipation;  Sea  spikes 


1.  Introduction 

Wave  breaking  has  profound  impacts  on  the  upper* 
ocean  turbulence  properties,  the  disruption  of  the  ocean 
surface  cool  skin,  the  generation  and  entrainment  of  air 
bubbles  into  the  water  column,  and  the  composition  of  the 
ocean  surface  roughness.  All  factors  cited  above  are 
known  to  produce  significant  modifications  of  the  gas 
transfer  rale  across  the  air  sea  interface  (c.g.,  Phillips. 
1085;  Wallace  and  Wirick,  1992;  Thorpe.  1993;  Melville. 
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1 994;  Monahan,  2002;  Asher ct  at,  2002).  Clarification  of 
the  spectral  properties  of  breaking  surface  w  aves  may 
contribute  to  a  better  understanding  of  the  gas  transfer 
process  in  the  ocean.  Wave  breaking  is  also  an  important 
subject  in  many  different  ocean  remote  sensing  applica¬ 
tions.  for  example,  sea  spikes  in  radar  scatter  from  the 
ocean  surface  (c.g.,  Lee  et  al..  1996;  Frasier  ct  aL  1 998; 
Liu  ct  al.,  1 998;  S’htlltps  el  al,  200 1 )  and  buhble  effects  on 
underwater  acoustics  {e.g.T  Ding  and  Farmer,  1994;  Deane 
and  Stokes,  2002;  Melville  and  Mutusov,  2002). 

In  this  paper,  experimental  results  on  the  length  scale 
of  breaking  waves  are  examined.  Section  2  describes  the 
indirect  approach  of  breaking  wave  analysis  through  the 
investigation  of  the  dissipation  function  of  the  wave 
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action  density  conservation  equation  (e.g.,  Phillips, 

1 984,  1 985),  The  method  has  been  applied  to  the  spectra 
of  intermediate-scale  waves  (from  0,02  to  6.3  m  long) 
collected  in  the  ocean  using  a  free-drifting  measurement 
technique  to  mitigate  the  problem  of  Doppler  frequency 
shift  in  converting  the  measured  encounter- frequency 
spectrum  to  the  wavenumber  domain  (Hwang  and 
Wang,  2004),  Interesting  features  of  the  spectral 
properties  of  intermediate-scale  surface  waves  in  the 
ocean  derived  from  the  analysis  arc  discussed.  In 
Section  3,  the  wind-speed  dependence  of  the  wave 
spectral  density  is  presented.  The  topic  is  of  great 
interest  in  many  areas  of  research  including  wave 
dynamics,  air  sea  interaction  and  microwave  remote 
sensing  of  the  ocean.  Results  from  earlier  studies  on 
surface  waves  of  length  scales  from  long  gravity  waves 
to  short  capillary  waves  are  found  to  be  in  good 
agreement  with  the  intermediate-scale  wave  data.  In 
Section  4,  the  dissipation  function  derived  from  the 
action  density  conservation  equation  is  described.  The 
most  interesting  result  is  that  the  dissipation  function 
displays  a  quasi-singular  property  in  the  middle  range 
(about  0,15  to  1,5  m)  of  the  intermediate-length  scale. 
Such  a  feature  is  inteq^reted  as  a  distinct  signature  of 
wave  breaking  in  wave  components  of  intermediate- 
length  scale.  Results  on  the  length  scale  of  wave 
breaking  obtained  from  more-direct  wave  breaking 
measurement  techniques,  including  tracking  of  bubble- 
generated  acoustic  noise  (Ding  and  Fanner,  1994),  radar 
sea  spikes  (Lee  ct  al.,  1996:  Liu  cl  a!.,  1998;  Frasier 
cl  al.,  1998;  Phillips  ct  al.,  2001),  whitecaps  (Melville 
and  Matusov,  2002),  and  high-speed  photographs  of 
bubble  plume  formation  (Deane  and  Stokes,  2002)  are 
investigated  in  Section  5.  It  is  shown  that  the  length 
scales  of  wave  breaking  derived  from  direct  and  indirect 
methods  are  in  very  good  agreement.  Additional 
discussions  on  issues  of  wave  dynamics  are  given  in 
Section  6  and  a  summary  is  presented  in  Section  7, 

2*  Analysis 

Wave  breaking  is  the  dominant  dissipation  mecha¬ 
nism  in  the  source  function  balance  of  the  action  or 
energy  conservation  equation  of  surface  waves.  The 
conservation  equation  of  wave  action,  A r(A),  can  be 
expressed  as  (e.g.,  Phillips,  1984,  1985) 


dJV 

d/ 


^  +  o. 

,  A 


(1) 


where  A  is  the  wavenumber  vector  with  modulus  A,  The 
three  terms  on  the  right-hand  side  of  (1)  represent  the 
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source  functions  due  to  nonlinear  wave  wave  interac¬ 
tion,  wind  input  and  breaking  dissipation.  N  is  related  to 
the  2D  wave  displacement  spectral  density,  x(A),  by  A' 
(A)=#x(A)/o\  g  the  gravitational  acceleration,  and  a  the 
intrinsic  frequency  of  the  wave  component.  In  the 
following,  the  result  applicable  to  the  omni-directional 
spectrum, %( A) %( k)k&Q  and  N( k ) _ gx( A)/a,  is 
discussed. 

Theoretical  and  experimental  studies  of  the  energy 
transfer  from  wind  to  waves  have  led  to  the  following 
parameterization  function  for  the  wind  input  source 
function  (Plant,  1982;  Phillips,  1984) 


where  m  =  0.04.  For  gravity  waves  several  times  shorter 
than  the  spectral  peak  component,  the  nonlinear  wave 
wave  interaction  term  is  much  smaller  than  the  other  two 
terms  (Phillips,  1977,  1984).  The  dynamic  balance  of 
intermediate- and  short-scale  gravity  waves,  therefore,  is 
mainly  determined  by  the  wind  input  and  breaking 
dissipation  terms.  Phillips  (1984,  1985)  suggested  the 
following  expression  for  the  omni-directional  dissipa¬ 
tion  sink  term 


D  —  gk~}f(B)<  (3) 

where  B  =  k^x  =  kJrrN/g  is  the  dimensionless  wave 
spectrum  representing  the  degree  of  saturation  (Phillips. 
1985).  Lq.  (2)  can  be  expressed  in  terms  of  B , 

=  gk  'm-  (4) 

At  equilibrium  dA7d/=0.  Equating  (3)  and  (4),  Phillips 
(1984)  concluded  that  the  solution  to  the  unknown 
dissipation  function  j\B)  relies  on  the  knowledge  of  the 
dependence  of  B  on  the  wind  speed.  Experimental  data 
show  that  the  dependence  of  the  spectral  density  on 
wind  speed  can  be  represented  by  a  power-law  function 
(Hwang  and  Wang.  2004) 


Assuming  a  power-law  relation  for  the  unknown 
function  in  the  dissipation  term  (Phillips,  1984,  1985) 


/(B)  -  AdBa* 


(6) 


the  expressions  for  c/d  and  become  simply 
tfd  =  1  +  — . 


(7a.  b) 
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3.  Spectral  dependence  an  wind  speed 

The  approach  described  above  is  most  suitable  for 
applications  in  the  region  of  intermediate-  and  short- 
scale  waves  where  the  magnitude  of  the  nonlinear 
wave- wave  interaction  term  is  much  smaller  than  that 
of  the  wind  input  or  the  breaking  dissipation  term* 
Obtaining  reliable  measurements  of  intermediate-  and 
short-scale  waves  in  the  ocean  is  a  complicated  task. 
One  of  the  most  difficult  issues  is  the  Doppler  frequency 
shift  that  renders  the  interpretation  of  the  length  scale  of 
measured  encounter  frequency  spectrum  an  uncertain 
task  (Hwang,  2006)*  This  problem  can  be  mitigated  by 
free-drifting  operation  that  provides  measurements  in  a 
frame  moving  with  the  adverting  currents  that  caused 
the  Doppler  frequency  shift.  Phillips  (1984)  method  of 
analysis  described  in  the  last  section  is  applied  to  the 
spectra  of  intermediate-scale  waves  (wavelengths  be¬ 
tween  0*02  and  6*3  m)  measured  in  the  ocean  using 
wave  gauges  mounted  on  a  free-drifting  platform 
(Hwang  and  Wang,  2004)*  The  results  show  that  B(u+.f 
c)  can  be  represented  by  a  power-law  function  (5).  An 
example  is  shown  in  Fig.  la* 

Recently,  laboratory  experiments  were  conducted  in 
the  large  wind  wave  facility  at  the  Institut  de  Recherche 
sur  les  Phenomenes  Hors  iiquilibre  (IRPHE),  Marseille* 
France.  The  wind-generated  surface  waves  are  measured 


(a) 


by  two  separate  sensor  systems.  The  first  system  is  a  two- 
dimensional  scanning  laser  slope  sensor  that  generates 
directly  the  wavenumber  spectrum  from  the  spatial 
images  of  the  surface  slopes  over  an  area  0.1  mx0. 1  m. 
The  second  system  constitutes  of  thin- wire  wave  gauges 
at  fixed  stations  that  provide  the  frequency  spectrum  of 
surface  waves  measured  at  essentially  the  same  wind  fetch 
as  that  of  the  scanning  slope  sensor  system.  The  labo¬ 
ratory  measurements  further  confirm  the  robust  power- 
law  relationship  of  B(ujc)  as  illustrated  by  the  examples 
shown  in  Fig.  lb.  From  these  field  and  laboratory  data,  Au 
and  oq  can  be  obtained  by  least -square  fitting  of  the 
spectral  measurements  for  each  wavenumber  component. 
The  results  are  given  in  Fig.  Ic  and  d,  respectively.  An 
interesting  feature  in  the  results  of  A o  and  ao  is  their 
nonmonotonic  dependence  on  k%  as  clearly  illustrated  in 
the  field  data.  The  results  from  laboratory  experiment  (all 
wind-sea  conditions)  show  a  similar  dependence  ot\4u 
and  eta  on  A-  in  the  wavenumber  range  comparable  to  those 
of  the  field  data,  but  because  the  length  scale  of  laboratory' 
wind-generated  w  aves  is  much  shorter  than  that  in  the 
field,  the  interesting  nonmonotonic  feature  revealed  by 
the  field  data  cannot  be  detected  in  the  laboratory  results. 
Also,  there  is  an  apparent  wavenumber  downshift  in  the 
c  luster  of  laboratory  data  derived  from  temporal  measure¬ 
ments  by  stationary  season;  when  compared  to  the  result 
derived  from  spatial  measurements.  This  is  likely  due  to 


(c) 


J  ig.  I .  Fxmuplcsofthc  power-law  dependence  of  BUt0k  y.  (a)  field  measurements  by  wave  gauges  mounted  on  a  free-drifting  platform,  (b>  laboratory 
measurements  by  stationary  wave  gauges  (temporal)  and  scanning  laser  slope  sensor  (spatial)  The  wavenumber  dependence  offe)  the  proportionality 
coefficient.  and  (d)  the  exponent,  ath  of  the  power-law  function  derked  from  both  field  and  laboratory  experiments 
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the  Doppler  frequency  shift  induced  by  the  orbital 
velocity  of  longer  waves.  Because  the  Jacobian  connect¬ 
ing  the  wavenumber-to-frequeney  conversion  is  nonlinear 
with  respect  to  the  adverting  current,  the  Doppler  fre¬ 
quency  shift  due  to  the  background  orbital  velocity  docs 
not  vanish  on  averaging  over  integral  cycles  of  long 
waves.  Instead,  the  advert  ion  by  orbital  velocity  produces 
a  net  downshift  in  the  resulting  encounter  frequency  for  a 
given  wavenumber  component  (Hwang.  2006). 

The  nonmonotonic  behavior  of  the  spectral  depen¬ 
dence  on  wind  speed  as  quantified  by  is  noteworthy. 
Discussions  of  this  feature  have  been  given  hv  Hwang  and 
Wang  (2004)  and  they  are  briefly  summarized  in  the 
following.  As  illustrated  in  Fig.  Id,  for  both  conditions  of 
wind  sea  and  mixed  sea,  the  wind  speed  dependence  of  the 
wave  spectral  density  approaches  linear  toward  the  long 
gravity  wave  region.  This  is  in  agreement  with  many 
earlier  investigations  on  the  wind-speed  dependence  of 
the  surface  wave  spectrum  in  the  equilibrium  range, 

yjk)  -  hu>g  =  A  \  (8) 

The  corresponding  dimensionless  spectrum  is 

BAk)  =  h(^j.  (9) 

The  range  of  b  varies  approximately  within  a  &ctor*of- 
two  range  based  on  field  observations,  as  lias  been 


summarized  by  many  researchers  (e.g„  Toba,  1973; 
Phillips,  1985;  ForristalL  1981;  Donctan  et  at.  1985; 
Hwang  cl  al.,  2000),  Based  on  the  2D  wavenumber 
spectrum  derived  from  the  spatial  3D  surface  wave 
topography  of  a  steady  and  equilibrium  wind-generated 
wave  field  obtained  by  an  airborne  scanning  lidar  system, 
Hwang  ct  al,  (2000)  reported  5.26x10  2.  The 
analysis  of  intermediate-scale  wave  spectra  also  yields 
similar  values  of  the  proportionality  coefficient  and 
exponent  (Aty  and  u(l)  in  the  lower  limit  of  the  resolved 
wavenumber  range  ( l  <  k  <  3 16  rad/ni)  ( Fig,  I  e-d). 

Banner  ct  al,  (1989)  reported  wavenumber  spectra  of 
surface  waves  derived  from  3D  surface  topography 
measured  by  stereo  photography.  The  resolved  wave¬ 
lengths  are  between  0.2  and  1 .6  m  (3 1  >  4  rad/m).  The 

range  of  wind  speeds  is  between  5.5  and  13*3  m/s.  Their 
results  show  a  very  weak  wind-speed  dependence  of  the 
wave  spectrum,  with  fl(A)  — h#*’  **  (indicated  by  a  short 
line-segment  in  Fig.  2).  The  analysis  of  intcrmediatc-scalc 
wave  spectra  expands  the  resolved  wavelengths  to  a  range 
from  0,02  to  6,3  m  ( 1  <  k  <3 1 6  rad/m).  For  wind  sea,  the 
wind- speed  exponent  drops  sharply  from  approaching  1 .0 
( linear)  for  long  gravity  waves  to  about  0.22  near  k --5  rad/ 
m;  th\  remains  to  be  less  than  0,5  through  A=4Q  rad/m.  In 
mixed  sea,  the  minimum  value  of  a0  is  also  about  0.22, 
and  fin  <0.3  for  k  between  8  and  25  rad/m,  and  *in  <0.5 
for  Ar  between  3  and  80  rad/m  (Fig,  2). 

For  short  gravity  w-aves  and  capillary  waves.  Held  data 
are  usually  obtained  by  radar  backscatter  measurements 


k  ( rjJ/m  I 

Fig,  2.  The  wind-speed  exponent  of  ii  given  wave  spectral  component  obtained  by  different  methods  including  radar  return  (Jones  and  Schmcder. 
197H;  Stewart,  1985;  Masuko  et  al  .  198fi;  Phillips,  1988;  Wcissman  cl  al.,  1 994;  Colton  el  al.,  1995:  summarized  in  Fig.  5  of  llwang.  1997),  stereo 
photography  (Banner  et  al..  1989),  and  free-dri  fling  nave  gauge  (l  twang  and  Wang,  2004);  kf  is  the  radar  wavelength 
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using  different  radar  frequencies  and  incident  angles  (c.g., 
Jones  and  Schroedcr,  1978;  Stewart,  1985;  Masukoct  a)., 
1986;  Phillips,  1988;  Weissman  ct  al.,  1 994;  Colton  et  a l, 
1995).  Fig.  5  of  Hwang  ( 1 997)  summarized  the  radar  data 
with  incidence  angles  between  25c  and  60°,  in  which 
range  the  dominant  scattering  mechanism  is  the  Bragg 
resonance  so  the  radar  cross  section  is  proportional  to  the 
spectral  density  of  the  resonant  roughness  spectral 
component.  The  radar  data  am  shown  in  Fig.  2  together 
with  the  result  derived  from  the  analysis  of  intermediate* 
scale  wave  spectra  described  above.  Despite  the  large 
scatter  in  the  radar  data,  there  are  apparent  similarities  in 
the  results  obtained  by  these  fundamentally  different 
approaches  of  obtaining  the  wind-speed  exponent  of  the 
wave  spectral  component, 

4,  Dissipation  function 

Fig.  3  shows  the  coefficients  (7)  of  die  dissipation 
function  (6)  obtained  from  field  and  laboratory 
measurements.  The  breaking  dissipation  function  dis¬ 
plays  a  quasi-singular  behavior  in  the  neighborhood  of 
meter-long  wave  components.  In  mixed  sea,  the 
breaking  scale  shifts  toward  smaller  breaker  size  and 
the  size  range  expands  in  comparison  with  those  of  the 
wind  sea,  reflecting  the  modification  by  long  waves  and 
interaction  between  background  orbital  velocity  and 


short  waves.  The  results  suggest  that  the  spectral 
signature  of  wave  breaking  is  localized  in  the  wave- 
number  space  in  the  intermediate-length  scale.  [Detailed 
breaking  wave  observations  also  show  the  obvious 
signature  at  the  dominant-wave  length  scale  (e.g,. 
Banner  et  al.,  2000,  2002).  The  results  presented  here 
do  not  resolve  wavelengths  near  the  spectral  peak 
component.  The  focus  of  the  present  research  is  on 
surface  waves  in  the  intermediate-length  scale.  These 
intermediate-scale  waves  contribute  more  than  78%  to 
the  total  mean-square  slope  of  gravity  waves  generated 
by  winds  up  to  20  m/s  (Hwang,  2005).]  From  the  poini- 
of-view  of  the  spectral  response  to  wind  input,  as  shown 
in  Fig.  Id  or  Fig.  2,  the  presence  of  a  localized  region  in 
the  wavenumber  space  where  the  spectral  density  is  only 
weakly  dependent  on  the  forcing  wind  condition 
suggests  a  localized  region  where  the  wave  growth  is 
quenched  by  strong  dissipation.  Further  discussions  of 
this  point  will  be  described  in  Section  6,  The  range  of 
wavelengths  with  oo<0.3  is  about  0.15  <A<  1.5  m. 
These  wave  components  arc  sufficiently  long  such  that 
Lhc  energy  loss  due  to  viscous  dissipation  (Lamb,  1945) 
or  generation  of  parasitic  capillary  waves  (Longuct- 
Higgins,  1992,  1995)  may  be  neglected.  It  is  reasonable 
to  assume  that  wave  breaking  is  the  major  dissipation 
mechanism  for  the  observed  peculiarity  in  the  dissipa¬ 
tion  function  of  intermediate-scale  waves. 


la)  (b) 


Fig,  V  The  wavenumber  dependence  of  (a)  ihc  proportionality  coefficient,  Aa-  and  (b)  die  exponent,  i/j.  of  the  power-taw  representation  of  the 
dissipation  function*  nil)  Results  from  free-drilling  measurements  m  the  field  in  wind  sea  and  mixed  sea  as  well  as  laboratory  data  by  stationary' 
temporal  and  spatial  measurements  are  presented. 
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Microwave  observations  of  the  ocean  surface  over 
the  Iasi  few  decades  indicate  that  Doppler  shifts  at  HH 
(horizontal  transmitting,  horizontal  receiving)  polariza¬ 
tion  and  low  grazing  angles  are  due  to  scatterers 
ira ve ling  at  the  speed  of  surface  waves  of  a  few  meters 
long  (e.g.,  Lee  el  al.*  1996;  Frasier  et  ah,  1998;  and  the 
references  therein).  This  feature  is  consistent  with  the 
strong  signature  of  wave  breaking  in  the  mtermediatc- 
length  wave  components.  In  the  next  section,  some 
recent  breaking- wave  observations  made  by  more-direct 
measurement  techniques  arc  summarized  to  compare 
with  the  result  using  the  indirect  approach  of  dissipa¬ 
tion-function  analysis. 

5*  Urea  king  length  scale 

In  the  last  decade  or  so*  many  more-direct  measure¬ 
ments  of  breaking  waves  have  heen  reported,  using 
techniques  such  as  passive  acoustic  tracking  of  the 
ambient  noise  produced  by  bubbles  generated  by 
breaking  waves  (Ding  and  Fanner.  1994),  radar  tracking 
of  sea  spikes  (Lee  et  al.,  1996:  Frasier  et  al.f  1998;  Liu 
ct  al.*  1998;  Pit i Hips  et  ah,  2001),  and  video  tracking  of 
whitecap  evolution  (Melville  and  Matusov,  2002). 
These  remote  sensing  techniques  usually  process  a 
very  large  population  of  breaking  events*  on  the  order  of 
tens  of  thousands.  From  these  analyses,  the  authors 
report  statistics  of  the  frequency  of  occurrences  of 


breaking  events  as  a  function  of  the  breaker  phase  speed, 
The  length  scales  of  breaking  waves  can  be  calculated 
from  the  reported  statistics,  they  range  from  about  0*1  to 
3  m  al  4  m/s  wind  speed  and  from  3  to  20  m  at  1 5  m/s 
wind  speed  (Fig.  4a).  Two  major  factors  contribute  to 
the  observed  wide  range  of  the  breaking  length  scales  at 
a  given  wind  speed.  The  first  is  the  dynamic  range  of  the 
sensors.  In  particular,  the  passive  acoustic  tracking  of 
breaking  events  cannot  detect  small-scale  breakings  due 
to  the  problem  of  ambient  noise  (Ding  and  Fanner, 
1994).  The  second  factor  is  the  stage  of  wave 
development.  The  experiments  reported  by  Lee  ct  al. 
(1996)  were  conducted  in  a  lake  and  a  protected  bay 
with  limited  wind  fetch  so  the  wave  field  is  relatively 
young  compared  to  the  wave  conditions  of  the  others 
obtained  in  the  open  ocean  environment.  The  second 
factor  can  be  compensated  somewhat  by  normalizing 
the  breaking  length  scale,  by  the  dominant 
wavelength  of  the  wave  field*  With  this  normaliza¬ 
tion,  the  majority  of  radar  data  are  clustered  in  the 
range  between  0*04  and  0.2  over  the  full  range  of  wind 
speeds  encountered  in  the  experiments  (Fig.  4b).  The 
mean  value  with  a  standard  deviation  is  X\Jkp- 0.092 ± 
0.051,  thus  the  representative  breaking  length  scale  is 
about  one  order  of  magnitude  shorter  than  the  dominant 
wavelength  of  the  wave  field  based  on  sea-spike 
observations.  [Although  the  normalized  breaking  length 
scale  seems  to  bring  together  data  sets  from  very 
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J  tgr  4.  (a)  The  representative  breaking  length  scale*  Xhm  obtained  by  different  more-direct  measurement  techniques  of  wave  breaking  including 
tracking  of  acoustic  noise,  whitecap*  and  sea  spikes.  Estimation  of  the  upper  and  lower  bounds  of  the  breaking  length  scales  from  the  dissipation 
analysis  is  graphed  as  boxes.  \  b)  Dimensionless  breaking  length  scale*  <v/T. 
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different  wave  development  stages*  it  should  be  pointed 
out  that  for  open  ocean  measurements  the  data  scatter 
with  breaking  length  given  in  dimensional  form  is  much 
smaller  than  that  given  in  dimensionless  form  (Fig*  4), 
The  cause  and  significance  of  ill  is  interesting  result  is 
not  obvious  at  this  stage  but  it  seems  to  suggest  that  the 
peak  wavelength  is  not  necessarily  the  proper  sealing 
factor  for  intermediate- length  wave  components  several 
times  shorter  than  the  spectral  peak  wavelength.  If  the 
scaling  length  is  the  wavelength  at  minimum  phase 
speed  (determined  by  the  ratio  of  the  gravitational 


acceleration  and  the  surface  tension  divided  by  the  water 
density*  and  is  almost  constant  in  sea  water  without 
surface  slicks)*  the  reduced  scatter  in  the  dimensional 
representation  of  the  breaking  length  scale  would  be 
sensible.  However*  as  discussed  earlier*  the  range  of  the 
length  scales  of  nonmonotonic  behavior  is  between  0.15 
and  1 .5  m*  and  tar  away  from  the  capillary  influence*  A 
better  explanation  remains  elusive*] 

The  breaking  length  scale  can  be  estimated  from  the 
indirect  breaking  analysis  using  the  dissipation  function 
by  setting  a  threshold  in  the  result  of  or  a*|(£) 


I  ig  5  Piuluics  from  high-speed  photography  of  bubble  plume  formation  under  breaking  waves  (|-ig.  2  of  Deane  and  Siokcs,  2002k 
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shown  in  Fig.  3.  For  example,  applying  a  threshold  of 
10*  lirncs  less  than  the  peak  value  oMti,  one  can  derive 
the  breaking  length  scales  from  the  lower  and  upper 
bounds  of  the  wave  numbers  that  A d  exceeds  the 
threshold.  The  result  is  shown  in  Fig.  4.  The  breaking 
length,  Ah,  is  between  O.h  and  1.5  m  for  wind  sea.  and 
between  0.15  and  1.2  m  for  mixed  sea.  The  ratio  X\JXP 
ranges  between  0.013  and  0.025  for  wind  sea,  and  From 
0.008  to  0.05  for  mixed  sea.  Similar  results  ean  be 
derived  by  using  a  threshold  (of  about  7)  in  f/j.  These 
values  of  breaking  length  scales  are  smaller  than  those 
derived  from  sea-spike  measurements.  It  is  likely  that 
sea-spike  tracking  overestimates  the  breaking  length 
scale  because  the  tracking  procedure  requires  a  finite 
dwell  time  of  the  breaking  event  in  the  field  of 
measurements.  For  example,  Frasier  et  at.  (1998)  set  a 
threshold  of  Us  minimal  dwell  time  for  the  sea-spike 
event  to  be  included  in  the  breaking  analysis,  thus 
excluding  shorter  breaking  events  in  the  final  statistics. 
This  problem  is  similar  to  the  ambient-noise  limitation 
on  the  length-scale  resolution  in  the  acoustic  tracking  of 
breaking-induced  acoustic  signature  encountered  by 
Ding  and  Farmer  (1994).  Taking  this  bias  into 
consideration,  it  is  concluded  that  the  results  from  direct 
and  indirect  measurements  of  the  breaking  length  scales 
are  comparable. 

Deane  and  Stokes  (2002)  conducted  laboratory  study 
on  the  generation  mechanism  of  bubble  plumes  under 
breaking  waves.  Using  high-speed  photography,  they 
captured  many  interesting  hydrodynamic  processes  in  the 
active  phase  of  wave  breaking.  Relevant  to  the  present 
study  is  the  sequence  of  pictures  depicting  the  develop¬ 
ment  of  bubble  plume  in  the  first  two  seconds  after 
formation  of  the  air  cavity  trapped  by  the  breaking  jet 
(their  Fig.  2,  reproduced  as  Fig.  5  here).  The  dominant 
frequency  of  the  wave  field  is  0.73  Hz  and  the 
corresponding  wavelength  is  2.3  m.  As  illustrated  by 
the  sequence  of  pictures,  after  more  than  one  wave  period 
following  the  formation  of  the  air  cavity  produced  by  the 
jet  of  wave  breaking,  the  length  scale  of  the  active 
breaking  region  with  intensive  dissipation  (the  vortex 
regions  in  the  photographs)  remains  on  the  order  of  0.1  m. 
This  gi\  es  an  estimate  of  the  ratio  X\ JXp  -0.025,  which  is 
in  very  good  agreement  with  the  results  derived  from  the 
analysis  of  the  dissipation  function  shown  in  Fig.  4b. 

6.  Discussions 

Investigation  of  the  breaking  dissipation  function  is 
of  great  interest  to  many  areas  of  research  ranging  from 
ocean  and  coastal  engineering,  wave  dynamics,  air-sea 
interaction,  remote  sensing,  and  numerical  wave 


models.  In  most  applications,  the  action  density 
conservation  equation  formulated  as  Eq.  ( l)  is  sufficient. 
In  coastal  areas,  the  bottom  dissipation  function  needs  to 
be  included  on  the  right-hand  side  of  the  equation.  For 
intermediate-scale  waves,  the  situation  is  not  as  clear.  A 
key  question  is  whether  wind  forcing  is  the  sole  source 
of  wave  generation.  The  results  from  radar  sea-spike 
analyses  are  very  instructive.  Many  of  ihose  experi¬ 
ments  arc  accompanied  with  detailed  video  recording  of 
the  wave  field  in  the  radar  field-of-vicw  (e.g.,  Frasier 
cl  al.,  1998;  Liu  et  al.,  1998).  The  close  association  of 
sea  spikes  with  (large-scale,  or  dominant)  wave 
breaking  or  steep  waves  is  convincingly  established, 
yet  the  phase  velocity  of  the  breakers  derived  from 
tracking  sea  spikes  is  much  slower  than  the  phase  speed 
of  the  dominant  wave  component.  The  logical  explana¬ 
tion  of  this  result  is  in  the  distinction  of  wave  breaking 
events  (as  in  dominant  wave  breaking)  vs.  breaking 
patches  (as  in  sea  spikes  or  the  length  scale  of  the 
dissipation  function).  The  breaking  patches  are  obvi¬ 
ously  generated  by  breaking  of  dominant  waves.  For  a 
very  short  moment,  they  retain  a  phase  speed  close  to  or 
exceeding  that  of  the  dominant  waves  (as  a  bound-wave 
component).  Judging  from  the  sea-spike  tracking 
analysis,  this  duration  is  very  short  and  for  the 
remaining  part  of  the  sea-spike  (breaking)  lifetime,  the 
breaking  patches  propagate  as  free  waves  with  phase 
velocities  according  to  their  sizes  (in  the  intermediate- 
length  scale).  In  other  words,  breaking  plays  double 
roles  in  the  intermediate-scale  wave  components:  it  is  a 
dissipation  function  of  wave  energy  and  ail  important 
source  function  of  imcimcdiatc-scalc  wave  generation. 

Although  it  is  reasonable  to  state  that  there  is  an 
additional  source  function  due  to  breaking-wave 
generation  in  the  intermediate-scale  spectral  compo¬ 
nents,  it  remains  a  difficult  task  explaining  the  lack  of 
response  to  wind- forcing  in  the  mid-range  components 
of  intermediate-scale  waves,  observed  by  both  in  situ 
wave  measurements  and  remote-sensing  radar  scatter 
analyses  (Fig.  2).  It  is  generally  accepted  that  breaking 
probability  increases  with  the  cube  of  wind  speed,  so 
incorporation  of  breaking  as  a  generation  term  for 
intcrmcdiatc-scalc  waves  should  increase  the  wind- 
speed  exponent  for  the  wave  components  affected, 
contradicting  to  the  observations  of  decreased  wind- 
speed  exponent  in  the  middle  range  of  the  intermediate- 
length  wave  components.  It  remains  a  tough  challenge 
to  explain  the  observed  peculiar  nonmonotonic  behavior 
of  the  spectral  properties  of  intermediate-scale  wave 
components  shown  in  Fig.  2. 

Eq.  (1)  is  obviously  inadequate  for  describing  the 
dynamics  of  intcrmcdiatc-scalc  waves.  While  it  worked 
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nicely  for  identifying  the  breaking  length  scales,  it  is 
emphasized  that  the  proposed  local  balance  of  breaking 
dissipation  and  wind  generation  functions  is  a  hypoth¬ 
esis  that  leads  to  die  prediction  of  enhanced  dissipation 
in  the  intermediate-length  wave  components.  Further¬ 
more,  in  side-by-side  comparison  of  radar  images  of  sea 
spikes  and  video  images  of  ocean  surface  waves,  it  was 
found  that  the  majority  of  sea -spike  events  arc 
associated  with  steep  waves  without  visible  breaking 
whitecaps  approximately  60%  tor  the  young  and  the 
developed  sea  and  92%  for  decaying  sea  (Liu  ct  a  I., 
1998).  It  is  not  clear  whether  ihcse  sca-sp ike-associated 
steep  waves  were  undergoing  microscale  breaking,  so 
an  unequivocal  causal  link  between  sea  spikes  and 
breaking  waves  is  yet  to  be  established.  The  observa¬ 
tions  described  above  arc  supportive  but  not  conclusive 
and  further  observations  of  enhanced  breaking  at  the 
intermediate-  and  short-scale  waves  arc  needed  to 
substantiate  this  prediction  conclusively. 

7*  Nummary 

Properties  of  wave  breaking  can  be  investigated 
through  the  dissipation  function  of  the  wave  action 
conservation  equation.  Following  an  approach  de¬ 
scribed  by  Phillips  (1984),  the  functional  dependence 
of  i?(ii*/c)  is  found  to  follow  closely  the  power-law 
relationship  based  on  data  obtained  from  both  Held  and 
laboratory  environments.  The  exponent  of  the  power- 
law,  represents  the  wind-speed  dependence  of  the 
spectral  wave  components.  Results  from  field  data  that 
resolve  wavelengths  between  0,02  and  6.3  m  reveal  a 
nonmonotonic  behavior  of  a p  as  a  function  of  k.  The 
spectral  densities  of  wave  components  between  0. 1 5  and 
1.5  m  long  depend  only  weakly  on  the  forcing  wind 
condition.  The  dissipation  function  displays  a  quasi- 
singular  behavior  in  the  corresponding  wavelength 
range,  suggesting  a  localized  region  in  the  wavenumber 
domain  with  a  strong  breaking  signature.  This  result  is 
consistent  with  extensive  observations  of  radar  sea 
spikes  during  the  last  few  decades  showing  a  close 
correlation  between  sea  spikes  and  ocean  surface 
scatterers  traveling  at  the  speed  of  surface  waves  a  few 
meters  long.  The  strong  signature  of  breaking  waves  in 
the  decimeter-  to  meter-long  waves  reflects  the 
important  impact  of  wave  breaking  on  the  dynamics  of 
intermediate-scale  ocean  surface  waves  that  are  the 
major  contributor  to  the  ocean  surface  roughness 
(Hwang,  2005).  The  close  correlation  between  the  gas 
transfer  velocity  and  the  ocean  surface  roughness 
(mean-square  slope)  has  been  demonstrated  repeatedly. 
A  better  understanding  of  the  properties  of  wave 


breaking  in  intcnncdiate-scalc  spectral  components 
can  lead  to  a  better  understanding  of  the  gas  transfer 
processes. 
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